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The nonintegrable Hamiltonian dynamics of particles placed in a symmetric, spatially periodic po- 
tential and subjected to a periodically varying field is explored. Such systems can exhibit a rich 
diversity of unusual transport features. In particular, depending on the setting of the initial phase 
of the drive, the possibility of a giant transient directed transport in a symmetric, space-periodic po- 
tential when driven with an adiabatically varying field arises. Here, we study the escape scenario and 
corresponding mean escape times of particles from a trapping region with the subsequent generation 
of a transient directed flow of an ensemble of particles. It is shown that for adiabatically slow in- 
clination modulations the unidirectional flow proceeds over giant distances. The direction of escape 
and, hence, of the flow is entirely governed whether the periodic force, modulating the inclination 
of the potential, starts out initially positive or negative. In the phase space, this transient directed 
flow is associated with a long-lasting motion taking place within ballistic channels contained in the 
non-uniform chaotic layer. We demonstrate that for adiabatic modulations all escaping particles 
move ballistically into the same direction, leading to a giant directed current. 
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I. INTRODUCTION 

Transport phenomena play a fundamental role in many 
physical systems. In this context, dissipative ratchet dy- 
namics has attracted considerable interest over the past 
years. Particles placed into periodic but sawtooth like 
potentials and driven by external forces or nonequilib- 
rium noise create a directed flow even if the forces and 
noise vanishes in average. Thus unbiased forces induce 
a directed motion, a concept which was successfully ap- 
plied to many different biological or mesoscopic systems. 
We refer here to the various overviews on molecular and 
Brownian motors Q-01- This physical concept was sub- 
sequently generalized to potential landscapes possessing 
reflection symmetry which in addition are subjected to 
asymmetrical driving [3|-[13| and, as well, to Hamiltonian 
transport in absence of dissipation and enduring agitat- 
ing fluctuations fill]-[27l|. 

The starting point for our investigation is a system 
obeying reflection symmetry, both in space and time. 
Being interested in the particle transport features in a 
nonintegrable situation we study the following equation 
of motion 

q = F (q) + F cos{Qt + 9 ) = F (q) + F(t) . (1) 

We consider spatially periodic forces Fo(q) = F (q + 1) 
such as e.g. the one provided by a completely symmet- 
ric potential Uo(q) — — cos(2irq) / (2ir) . The time-varying 
ac-force in |T|), F(t), periodically modulates the inclina- 
tion of the space-periodic potential Uo(q). This set-up 
destroys the integrability of the system dynamics. In 



particular, around the corresponding separatrix in phase 
space a chaotic layer develops. Like in multistable poten- 
tials the motion becomes irregular and trajectories jump 
erratically from one potential well to another, being not 
always the adjacent one. Hence particles are scattered 
by the nonlinear forces and, obviously, the property of a 
directed current will typically be lost. 

Remarkably, as pointed out in prior literature [la]. [la], 
in the system Jl]) there results an (unexpected) asymme- 
try of the flux of particles, emanating from one poten- 
tial well, and flowing to the left and right potential wells 
which indicates the existence of directed transport with- 
out breaking the reflection symmetry in space and time 
in this system. One reason for the occurrence of phase- 
dependent directed transport is the lowering of the sym- 
metry of the flow in phase space by the ac-field where 
this asymmetry vanishes only for specific values of the 
initial phase ©0 [la] . 

Also of interest in this context is Ref. [28| . Therein the 
authors report further on this exceptional situation and 
show that directed transport is sustained on fairly long 
time scales despite the presence of chaos. In particular it 
has been demonstrated that for sufficiently small forcing 
frequencies, fl -C 1, the width of the arising chaotic layer 
diverges leading to a strong enhancement of the chaotic 
transport (28)] . 

Here we complement and extend those prior studies 
[la], [28j to the problem of escape and successively main- 
taining exclusive directed motion. More precisely, parti- 
cles from a large ensemble that are initially contained in a 
well of the potential have to surmount the corresponding 



barrier. We show that adiabatically slow modulations of 
the potential lead to that all of the escaping particles not 
only leave the potential well in the same direction but 
subsequently enter the regime of long-lasting transients 
for which the transport proceeds unidirectional in a bal- 
listic fashion. The direction of this arising flow over giant 
distances depends on the initial phase 60 about which we 
do not average. 

That trajectories contained initially in the interior of 
the separatrix can escape to neighboring wells of the po- 
tential is related with sweeping across the chaotic layer 
and crossing the separatrix. Precisely for this situation 
the authors in []j| deomstrated (applying ac-forces with 
frequencies fl ~ 0(1)) that the mean time-averaged ve- 
locity of a particle ensemble decays inversely proportional 
with time and thus tending to zero asymptotically. In 
contrast, we demonstrate here that for driving with a 
sufficiently slow ac-field the mean time-averaged veloc- 
ity attains a finite value, virtually not altering on long 
time scales. As we shall show this is due to the fact 
that there occurs only a single event of separatrix cross- 
ing, namely the one guiding a particle from the inner to 
the outer region. Despite that afterwards the particles 
closely re-approach the separatrix periodically in time 
further crossings are avoided. 

We further provide reasonings for the occurrence of the 
flow on the basis of the underlying phase space structure. 
The only assumption we have to make is that particles 
which are initially trapped in the interior of the separa- 
trix are included in the arising chaotic layer around the 
separatrix in phase space. We will find and discuss the 
situation that large ensembles of particles do not only 
escape from the separatrix but also move in the same 
direction. This in turn yields a giant current. 

The paper is organized as follows: In the next sec- 
tion we introduce the model of the particle motion in a 
symmetric, periodic potential, associated with the force. 
The inclination of this potential is being time-periodically 
modulated. The structure of the phase space is eluci- 
dated. In section IIIII we investigate the influence of the 
modulation frequency on the formation of stochastic lay- 
ers in phase space. In particular, we determine the range 
in which the momentum variable can vary. The analysis 
in section IIVI concerns the chaos- induced escape of indi- 
vidual trajectories from the interior of the separatrix. In 
the subsequent section [V] an explanation for the effect of 
an enormous enhancement occurring in the adiabatically 
driven system is given in terms of the underlying phase 
space dynamics. The dynamics of ensembles of particles 
contained initially in the interior of the separatrix and 
their contributions to a directed flow are addressed in 
section IVII 

We close with a summary of our obtained results. 



II. THE FORCED NONLINEAR OSCILLATOR 
MODEL 

The equation (JTJ) , as a driven nonlinear Hamiltonian 
system in one dimension, is derived from the following 
Hamiltonian in (dimensionless) form 

2 2 

H = ^- + U(q,t) = ^ + U (q) + U 1 (q,t). (2) 

Therein p and q denote the canonically conjugate momen- 
tum and position of a particle evolving in the periodic, 
spatially-symmetric potential of unit period, i.e., 

U (q) = U (q + 1) = - cos(27rg)/(27r) . (3) 

An external, time-dependent forcing field 



Ui (q, t) = -F sm(il t + Q )q 



(4) 



serves for time-periodic modulations of the inclination 
of the potential Uq. We underline that the system is 
unbiased in the sense that the force averaged over time 
and space vanishes, i.e. 



T=27r/n 



dq 



dffiUo. 



(5) 



For a static inclination, i.e. for U\{q) = —Fq the sys- 
tem dynamics is integrable and the solutions are con- 
tained in the level set 



H 



P 



1 

2^ 



cos(27rq) - Fq = E" 



(6) 



There exist unstable saddles at q^ = 0.5 + k — 
arcsin(F)/(27r) and stable centers at q* = k — 
arcsin(i 7 ')/(27r) with integer values k = 0, ±1,±2 .... In 
the non-inclined case, F = 0, neighboring hyperbolic 
points are linked via heteroclinic connections and with 
inclination-modulation F^O each hyperbolic point is 
linked with itself via a homoclinic connection as illus- 
trated in Fig. [TJ According to the location of the saddles 
to the right (left) of the centers we denote the chain of 
homoclinic connections in the left (right) panel as right- 
oriented (left-oriented). For later reference, we point to 
the open channel arising between two neighboring separa- 
trix loops for inclined potentials. The separatrix energy 
attributed to the saddle at q^ is given by 



separatrix r\ 



y/l -F 2 - F[(0.5 + fc)27r - arcsin(F)] 

(7)" 
The maximal extension in momentum of the separatrix 

is 



pL = ± - Vi - f 2 

1 7T 



1/2 



(«) 



and the separatrix crosses the q— axis (turning point of 
the corresponding homoclinic orbit) at a value qo that is 
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Figure 1: (Color online) The separatrix in the phase plane. 
Top panel: A chain of homoclinic connections (solid line) for 
the potential with negative inclination due to a static force of 
amplitude F — 0.05. The dashed line represents the chain of 
heteroclinic connections belonging to a non-inclined potential. 
Bottom panel: Same as in the top panel except that now the 
inclination is reversed, i.e. F — —0.05, implying a positive 
inclination. The visible small openings in the solid lines of the 
separatrix refer to the open channels for ballistic transport. 



determined by the solution of the transcendental equa- 
tion 



F° 

separatrix 



1 

2^ 



cos(27rg ) + Fq = . 



(9) 



Note that the separatrix loops for F = completely 
comprises those for f ^ 0, For time-periodic modula- 
tions of the inclination which are imposed by the time- 
dependent potential U\{q,t), the phase space structure 
becomes governed by breathing of the separatrix loops 
where the left-oriented and right-oriented phase space 
structures displayed in Fig. Q] represent the "turning 
points" of the breathing. In-between these turnings the 
area enclosed by a separatrix loop is periodically chang- 
ing and the minimum and maximum is obtained when 
fit + O = 7r/4, 37r/4 and fit + O = 0, n, 2n respectively. 
Moreover, with the application of a time-dependent field 
U\ (q, t) a breaking of the integrability of the dynamics 



Figure 2: (Color online) Stroboscopic Poincare-plots for an 
angular frequency fl = 0.1, force strength F = 0.05 and initial 
phase So = 0. The top and the bottom panel are taken 
at tk = ( ^ + fc) T (maximal negative inclination) and tk = 
(| + k) T (maximal positive inclination), respectively. The 
coordinate q is represented mod 1. 



is expected. In particular, around the separatrix of the 
unperturbed system a chaotic layer forms. These oscil- 
lations of the separatrix between the left-oriented and 
right-oriented structures have to be distinguished from 
the pulsations of the width of the pendulum separatrix 
considered in [23, [30, uLU] ■ In Ref . [32] the chaos dynam- 
ics in a system with periodically disappearing separatrix 
was considered. 



III. CHAOTIC LAYERS AND 
POINCARE-PLOTS 

To illustrate the influence of the angular driving fre- 
quency on the behavior of the system we depict in Figs. El 
and [3] the stroboscopic Poincare-plots at successive pe- 
riods of the driving term T — 2tt/Q in the p — q— plane 
for force strength F = 0.05, phase 0o = and two dif- 
ferent driving frequencies fl — 0.1 and fl — 0.01. The 
simulation time is chosen such that fit = 10 4 holds and 
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Figure 3: (Color online) Stroboscopic Poincare-plots for an 
angular frequency fi = 0.01, force strength F — 0.05 and 
initial phase Go = 0. The top and the bottom panel are 
taken at tk — (^ + k)T (maximal negative inclination) and 
tk = (i + k)T (maximal positive inclination), respectively. 
The coordinate q is represented mod 1. 



therefore, the number of periodic changes 10 4 /ir of the 
inclination of the potential in the simulation time interval 
is the same for any frequency. 

The top panel in Figs. [2] and [3] corresponds to stro- 
boscopic plots at tk = (1/4 + k)T with k = 0,1,2,... 
where the potential assumes a maximal negative inclina- 
tion. The bottom panel belongs to stroboscopic plots at 
tk = (3/4 + k)T, i.e. the potential has a maximal posi- 
tive inclination. For fairly fast varying modulations with 
Q = 1 a chaotic layer of small width develops around 
the separatrix, whereas the motion remains trapped and 
regular within the large island of stability with its center 
at the origin of the phase plane. The deformed hori- 
zontal lines above and below the large island of stability 
correspond to KAM tori which act as barriers for trans- 
port impeding larger upwards and downwards excursions 
of the momentum variable. Note that the extension of 
the chaotic layer in momentum direction remains equal 
regardless of the sign of the inclination of the potential. 

Nevertheless, for slow modulation with fl — 0.01 this 



picture changes drastically. Many KAM tori become de- 
stroyed and the only surviving ones lie in the region of 
large \p\ > 5.8. There remains still an island of stability 
corresponding to bounded and regular motion inside the 
potential well. On the other hand, the chaotic layer has 
grown considerably in momentum direction compared to 
the previous case of il — 1 . Most strikingly, independent 
of the sign of the inclination of the potential the density 
of the points in the stroboscopic plots is much higher in 
the region of positive momentum than in the region of 
negative ones. Thus, the momentum variable is allowed 
to raise to fairly large positive values. Due to symmetry 
for an initial phase value Oo = 7r of the external field, 
i.e. when during the first half-period < t < T/2 the 
inclination becomes positive, an equivalent behavior is 
observed except that the cloud of points now penetrates 
more into the range of negative p. 

In the following we focus our interest on the generation 
of a directed flow for a large ensemble of particles where 
all of them perform exclusively ballistic motion in the 
same direction. Note that this represents a far stronger 
condition than merely observing directed, diffusive trans- 
port of particles as illustrated in [15j. Our set up is the 
following: If the potential assumes, let us say, a nega- 
tive inclination due to a static force —dUi/dq — F we 
suppose that the dynamics of the particles is bounded in 
one well and the corresponding potential barrier is insur- 
mountable, i.e. for all particles it holds that E° ticle — 
H = p 2 /2 + U {q) < E° separatrlx . Accordingly, the initial 
conditions for our simulations are distributed in the inte- 
rior of the corresponding separatrix loop (see also below 
in Fig. [8]). Obviously, with static inclination of the po- 
tential the integrable dynamics is characterized by oscil- 
lations around the stable elliptic center in the separatrix 
loop and hence, the particles remain trapped in the po- 
tential well. However, this scenario drastically changes 
when the time-dependent perturbations destroy the inte- 
grability. In particular homoclinic chaos is present in the 
driven dynamics. In fact, those trajectories seized by the 
arising chaotic layer may manage to escape from the inte- 
rior region of the broken separatrix. We emphasize that 
this trapping-detrapping transition can only be triggered 
by the chaotic dynamics within the chaotic layer. 

When applying the time-periodic modulation of the 
inclination with initial phase set at 8o = our starting 
point is a non-inclined potential. Motion is then sup- 
posed to proceed towards the right. Since in the initial 
stage < t < T/2, the inclination passes from zero value 
to its maximal negative value it contributes to a contin- 
ually growing positive force which promotes the desired 
motion towards the right. The question then is: What 
determines an efficient escape of trapped particles start- 
ing out from distributed initial conditions? Is the direc- 
tion of the motion of the escaped particles determined 
completely by the choice of the phase of the modulation 
term? If so, for which parameter values do the escaped 
particles keep moving (on average) in the preferred right 
direction on longer time scales ? 



IV. SEPARATRIX CROSSING AND ENERGY 
VARIATIONS 

In this section we discuss the trajectories of particles 
escaping from the interior of the separatrix loop. Fur- 
thermore, the possibility of energy growth after separa- 
trix crossing for slow modulations of the inclination is 
studied. 



Separatrix crossing and escape 

In Figs.[4]we depict the role of the angular frequency on 
the time evolution of trajectories for an initial condition 
contained in the region in the interior of the unperturbed 
separatrix in which the chaotic layer arises. For moder- 
ate and intermediate angular driving frequencies Cl = 1 
and £1 — 0.1, respectively, the coordinate q(t) behaves 
chaotically, exhibiting sudden and unpredictable changes 
of direction. 

In clear contrast, for slow driving at SI = 0.01 the 
coordinate dynamics not only seemingly behaves more 
regular but also grows on average upon evolving time. 
Due to symmetry the direction of the motion is reversed 
by changing the initial position q(0) to —q(0) and taking 
for the phase 9o = n. 

Furthermore, there are alternating long and short in- 
tervals during which the particle moves straightforwardly 
towards the right and left, respectively. Correspondingly, 
the momentum p(t) evolves in phases with positive value 
(motion towards the right) that are longer lasting than 
the phases when the momentum is negative (leftwards 
motion) (see also further below in Fig. \7$) . As Fig. [4] 
also reveals, the directed motion is maintained on a very 
long time scale where the coordinate assumes huge val- 
ues. Further details concerning the time scale of unidi- 
rectional motion are contained in Sec. IVI1 In the bottom 
panel of Fig. [4] one recognizes that for times t > 10 3 
the trajectory is trapped in a ballistic channel where the 
particle behaves effectively like a free particle propagat- 
ing ballistically. Notice that the simulation time interval 
T s = 2 x 10 6 is equivalent to almost 8 x 10 5 and 3184 
times the period duration for harmonic oscillations near 
the bottom of a potential well and the external modula- 
tion, respectively. 

Such motion in a ballistic channel occurs for the com- 
plex dynamics of systems with a mixed phase space [33j]- 
|3a |. In more detail the (broad) chaotic layer is not uni- 
form and contains cantori which can severely restrict the 
transport in phase space and thus effectively partition 
the chaotic layer bounded from below and above by non- 
contractible KAM tori [34J]-J36J]. Islands of regular motion 
that are situated at the upper and lower boundary of the 
layer in the vicinity of the confining KAM tori possess 
non-zero winding numbers and thus facilitate transport. 
The motion around these islands is characterized by the 
stickiness to them [3J|-[30|, that can lead to trapping of 
the trajectory for a long time resulting in ballistic mo 
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Figure 4: (Color online) Time evolution of the coordinate 
q(t) depicted with panels (a-c) for three different angular 
driving frequencies f2 (indicated in the corresponding pan- 
els (a)-(c)) but equally chosen initial conditions p(0) = and 
q(0) = 0.423. The remaining parameter values are F = 0.05 
and 0o = 0. In the central panel (b) we also depict the trajec- 
tory for initial coordinate q(0) — —0.423 and phase 8o = n for 
which the motion proceeds in the negative direction. In the 
bottom panel (c) the long-lasting ballistic motion of the par- 
ticle is illustrated. Note the large values for the simulation 
times > 10 6 in panel (c). The inset depicts a stroboscopic 
Poincare plot (at tk = Tk) of the trajectory, revealing the 
motion in a ballistic channel. 



tion [371. l38l|. For some islands the sticking times to the 
boundary of the islands can be anomalously long [39LI40I]. 



Energy variations 



It is illustrative to consider the difference, 



AE = E, 



particle 



between the particle energy 



E 



separatrix 



(10) 



EparUcle = y + U (q) - Fsin(0)g , (11) 

and the energy of the "frozen" separatrix at time t 

E S e P aratrix = U (q*) - F Sm(9)<^ , (12) 

where 9 = Vtt + 9 and gj = 0.5 + k + 
arcsin(i 7 'sin(0))/(27r) denote the instantaneous position 
of the corresponding hyperbolic point when the trajec- 
tory traverses the actual range 



Qt 1 <Q<Qs 



if p>0 

if P <0. (13) 



(Note that as q is a dynamical variable so is k.) 

As Fig. [5] reveals, for a fast modulation f2 = 1 the 
energetic difference AE changes frequently the sign, cor- 
responding to the trajectories' repeated leaving and re- 
entering of the interior of the instantaneous separatrix. 
Nevertheless, in-between separatrix crossings, the par- 
ticle can be trapped in a potential well for some time. 
Since AE stays close to zero the trajectory remains close 
to the separatrix for most of the time. Decreasing the 
angular driving frequency to ft — 0.1 has the effect that 
the number of separatrix crossings diminishes. Occasion- 
ally there appear interludes during which the trajectory 
escapes from the separatrix region giving rise to consid- 
erable coordinate changes q (cf. Fig. [4]). However, these 
changes are not coordinated, so that in essence no di- 
rected motion results. 

This behavior drastically differs from that occurring 
for a slow modulation Q = 0.01. Once the trajectory 
has crossed the separatrix it gathers momentum as long 
as the inclination remains non-positive. This goes along 
with an increase in the energetic difference AE. After- 
wards, during the depreciation period, i.e. when the in- 
clination of the potential is positive, the trajectory moves 
back towards the instantaneous separatrix, but never ac- 
tually re-enters its interior region. Notably, we followed 
the evolution on time scales as long as 10 7 where this 
behavior still holds true. In other words reversions of the 
direction of motion proceed in an open channel between 
two adjacent separatrix loops of the inclined potential as 
long as the adiabaticity condition 
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Figure 5: (Color online) Temporal behavior of the energy 
difference AE as defined in (| 10 j) for three different angular 
driving frequencies $1 (as indicated in the top, central and 
bottom panel) but equally chosen initial conditions p(0) = 
and q(0) — 0.423. The remaining parameter values are 
F — 0.05 and 6o = 0. In the inset in the bottom panel one 
recognizes that only a single separatrix crossing takes place. 



is satisfied, where for near-separatrix motion the period 
duration T tra jectory is asymptotically determined by 



- trajectory 



/2tt 



: 111 



AE» 



- E° I 

particle separatrix I 



(15) 

For the evolution depicted in the bottom panel in Fig. [4] 
and Fig. [5] this condition <fl4|) is obeyed. 
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Figure 6: (Color online) Illustration of the escape from the re- 
gion of the hyperbolic point in the phase plane. The particle 
gathers further momentum after the final libration in the in- 
terior of the separatrix has taken place (oscillatory solid line). 
The straight dashed lines run in the directions of the eigen- 
vectors belonging to the eigenvalues of the system linearized 
about the hyperbolic point. Sufficiently close to the hyper- 
bolic point they represent the stable and unstable mani fold 
W a and W u . The two horizontal solid lines at p± = ±^j2/ir 
confine the range of the momenta variations of the unper- 
turbed separatrix. The dashed closed line to the right of the 
hyperbolic point represents the instantaneous separatrix loop 
frozen at the moment when p attains the value p+. The pa- 
rameter values are given by F = 0.05, Q, — 0.01 and 6o = 0. 



The influence of the period of the driving force on the 
escape of the trajectory through the saddle point region 
can be elucidated as follows: Near the saddle point q s 
the dynamics is given by 



d*q 

dt 2 



+ U[l(q s )q = Fsin(nt-l-eo). 



where 



ugh) = 



d 2 U 
dq 2 



= —2tt — —a . 



\q=q 



(16) 



(17) 



The solution of equation (fT6| with initial condition p(0) 
and q(0) is given by 



P(t) 



p(0)+Fcos(@ 
q(0) + F 



fi 



a + n 2 

sin(9o 
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(18) 
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Figure 7: (Color online) Illustration of the enhanced motion 
at work for directed transport. Shown is the time evolu- 
tion of the momentum p and the corresponding force term 
F(t) = F sin(fit) with amplitude F = 0.05. For the sake 
of comparison F(t) is multiplied by a factor of 160. The re- 
maining parameter values are Q, = 0.01 and Oo = 0. Further 
details are discussed in the text. 



For directed motion to occur it is important that the 
trajectory, after having crossed the separatrix, gathers 
enough momentum that a sufficient distance to the sep- 
aratrix attributed to the saddle point (s) of the unsta- 
ble equilibrium of the next adjacent potential well(s) is 
reached before the inclination of the potential is reversed. 
Such behavior is illustrated with Fig. [6] for slow driving 
with 17 = 0.01. 

From the behavior of the solution ljT8")l . (fT9]) with ini- 
tial conditions situated close to the hyperbolic point we 
find that for relatively fast driving with ft = 1 the dis- 
tance of the trajectory immediately grows because of the 
still rapidly decreasing inclination of the potential being 
connected with increasing momentum and coordinate of 
the particle. Rather soon the sign of potential inclina- 
tion, and thus, the direction of the particle motion are 
reversed while the trajectory has not departed from the 
region close to the hyperbolic point. Correspondingly the 
trajectory, whilst being still in the neighborhood of the 
hyperbolic point, approaches the nearby part of the adja- 
cent separatrix (represented by the dashed line in Fig. [6]) 
and eventually crosses it soon after the escape. Contrar- 
ily, for slower modulations SI < 0.1 the trajectory slowly 
but continually increases its distance to the nearby part 
of the adjacent separatrix loop. In the vicinity of the hy- 
perbolic point the solutions (fl8|l . (fT9|) with fl < 0.1 reflect 
this behavior. 



V. MOMENTUM GROWTH 

Here we consider the mechanism for the speed-up of a 
particle that has escaped from the separatrix and derive 
an expression for the effective gain in momentum. 



Forced speed-up of particles 

To gain better insight into the mechanism that causes 
the speed-up of the particle motion and the transition 
from librations to rotations we depict in Fig. [7] the time 
evolution of the momentum for the driven dynamics to- 
gether with the corresponding force F(t) = —dUi/dq = 
Fsia(flt) with ft = 0.01. Up to times t < 132 the trajec- 
tory performs librations in the potential well. During the 
evolution the particle adiabatically increases its momen- 
tum so that with each turn in the potential it approaches 
closer the saddle point. Eventually, a separatrix cross- 
ing takes place when the trajectory is in the vicinity of 
the hyperbolic point where the motion becomes highly 
irregular, due to the intricate dynamics connected with 
the transversal intersections of the invariant manifolds 
|4ll . l42l |. Subsequently to such an escape, the motion 
proceeds along the direction of the unstable manifold of 
the hyperbolic equilibrium as further illustrated in Fig. El 
Most importantly, the separatrix crossing occurs at an 
instant of time t escape < T/4, for which the negative 
inclination has not yet reached its maximum. During 
the ongoing particle motion towards the right the mo- 
mentum variable performs small-amplitude oscillations, 
while its average value gr ows monotonically. Eventually, 
the largest value p = i/2/tt on the unperturbed sepa- 
ratrix (indicated by the upper straight line in Fig. El is 
exceeded. As the inclination of the potential assumes 
increasingly negative values the particles' momentum is 
raised further. At the end of the enhancement period, 
designated by Tf = T/2 — t escape in Fig. [J, the mo- 
mentum has grown to a maximum value of p m ax = 8. 
Subsequently, when for t > T/2 the inclination of the 
potential is positive, the particle is in the phase of depre- 
ciation, denoted by Tj.. Then, the momentum is reduced 
steadily and reverts to zero value at an instant of time 
tretum = T/2 + Xj > 3T/4. Due to symmetry, one ob- 
tains that Tf = Tj.. 

However, after reversal, i.e. when p < 0, there remains 
only comparatively little time, namely Ty, — T—t return < 
T/4, during which the motion to the left is enhanced. 
Consequently, the asymmetry in the enhancement and 
depreciation phases, viz. the fact that T/2 — t escape > 
T/4 > T — t re t U rn, serves for a rather long period of right- 
wards motion compared to the leftwards motion. There- 
fore the effective motion of the particle proceeds to the 
right. 



Momentum gain 

In fact the momentum gain for a particle exerted to the 
force F(t) = F sm(Qt + ©o) can be estimated as follows: 
Without loss of generality we consider phases in the in- 
terval < Go < 7r/2 leading to large motion towards pos- 
itive momenta if the particles escape at instants of time 
tescape such that Qt escape + O < 7r/2 is satisfied. Then 
there remains the time interval t escape < t < (ir — Oq)/U 



during which the particle still experiences an enhance- 
ment in the right direction. For times (ir — Qq)/Q, < t < 
(2tt — ©o)/^ the force acts in the opposite direction. In 
particular for t return < t < (2n — <do)/Q the momentum 
evolves with its sign reversed compared to the previous 
enhancement period. One obtains then 



Ap 



<.(7i--e )/fi / .(27r-e )/o 



drp 



= / + / 

\Jt C3capc Jt return ) 

x dr [- sin{2nq) + F sin(Or + 9 )] ■ (20) 

Due to symmetry, it holds that t return — (2ir — 6o)/^ — 
tescape- Furthermore, for small fl the oscillating part con- 
nected with the first term in the integral averages to zero 
on the time scale t escape < t < (2ir — Qq)/^ and we find 



F 
Ap = 2— cos(ilt escape + 6 ) . 



(21) 



In general, the smaller ft the higher is the gain in momen- 
tum (see also [28j]). In principle, for a sufficiently small 
frequency fl the gain can become arbitrarily large. Ap is 
non-negative only if 



te 



< 



f-e 



n 



(22) 



Integrating over the interval of escape times with the 
phase 8o being held fixed yields the integrated momen- 
tum gain 



< Ap > = 



(7r/2-e )/n 



dt. 



escape 



^P\J' escape ) 



= 2 — [1-sinen]. 



(23) 



In conclusion, the integral momentum gain is at its max- 
imum at 6o = and diminishes monotonically towards 
zero at Oo = 7r/2, underlining the vital (symmetry- 
breaking) role of the initial phase Oo for the enhancement 
process (see also |15l|). 



VI. ENSEMBLE DYNAMICS AND CURRENT 

In this section we investigate the behavior of the driv- 
ing induced current for an ensemble of trapped parti- 
cles. Their initial conditions are distributed such that 
if the potential had the static inclination F the associ- 
ated trapped trajectories cannot cross the corresponding 
separatrix loop. Likewise, as done in section IIVI for the 
chaos-induced detrapping, we apply a time-periodic mod- 
ulation of the inclination where the potential is initially 
non-inclined, i.e. we use a fixed initial phase ©o = 0. 
The main objective is to demonstrate that for imposed 
slow modulations a preferred direction of motion emerges 
although on average the potential landscape remains un- 
biased. 



a. 




Figure 8: (Color online) The separatrix loop of the non- 
inclined potential (dashed line) and for the inclined potential 
due to a static force set at F — 0.05 (solid line). The escape 
set comprising those initial conditions that lead to escape out 
from the solid separatrix loop is given by the hatched area (in 
green online). 
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Figure 9: (Color online) The relation between the particle 
energy E particle and the escape time for three different fre- 
quencies as indicated in the plot. The vertical dashed lines 
marked as (1),(2) and (3) are at the positions of the enhance- 
ment boundary T/4 = 7r/(2fi) for $7 = 1, marked by (1), 
ft — 0.1, indicated as (2), and fl — 0.01, marked as (3), re- 
spectively. The remaining parameter values are F — 0.05 and 
G o = 0. 



A. Escape times 

First, we identified the escape set, comprising of all 
those initial conditions that are contained in the interior 
of the static separatrix loop that lead to chaos-induced 
escape. To this end, the separatrix loop has been pop- 
ulated densely with points corresponding to initial con- 
ditions for the dynamics of the periodically driven sys- 
tem. The resulting escape set is displayed in Fig. [8] by 
the hatched area. We recall that escape, i.e. the de- 
trapping, requires a separatrix crossing. Apparently, the 
moment of the first separatrix crossing, i.e. the trapping- 



detrapping transition, governs the efficiency of the en- 
hancement process, as discussed in section ITVl Of great 
importance for the speed-up of an escaped particle is that 
the first separatrix crossing occurs at instants of time 
tescape < T/4 = 7r/(2fi), which determines the bound- 
ary for enhancement. More specifically, according to 
Eq. l[2"Tj) . the length of the enhancement period deter- 
mines the growth in momentum and thus the value of 
the resulting current. 

From Fig. [9] we deduce the corresponding time scale 
for escape beyond the separatrix. While for the largest 
angular driving frequency 0=1 the escape takes place 
at times far beyond the boundary, marked as (1) in Fig. 
[9j we note that for an intermediate angular frequency 
at O = 0.1 almost all escape events occur before the 
corresponding enhancement boundary T/4, marked as 
(2). For slow driving at $1 = 0.01, however, the scales 

Of tea 



. stretch over a wide range, but practically all 
escape events do occur at t escape < T/4, marked as (3). 
Moreover, due to the irregular nature of the underly- 
ing dynamics, for an ensemble of escaping particles with 
distributed initial conditions the moments of their corre- 
sponding first separatrix crossing depend sensitively on 
the initial conditions. Trajectories starting out very close 
to the hyperbolic point are the first to escape, whereas 
for those initial conditions which are located away from 
the separatrix considerable time passes until escape takes 
place. From Fig.[l0l depicting the mean escape time ver- 
sus the angular driving frequency Q, we infer that for 
low angular frequencies fl < 0.04 the escape times are 
sufficiently smaller than T/4. In contrast, for $1 > 0.04 
at least the upper error bars of the escape times are of 
the order of, or higher, than T/4. This hampers a pro- 
nounced enhancement and in turn tends to suppress the 
directed flow. Connecting the lower error bars of the 
mean escape time by a line we find that the latter is very 
well fitted by the expression T m i n = l/(2f2). Since this 
line is parallel to the one of the mean escape time the 
latter obeys the l/f2-dependence as well. 



B. Averaged directed flow 

With regard to directed flow we note that for $1 < 0.04 
we find that for the entire escape set the associated tra- 
jectories move in the same direction over an extended 
period of time. Thus, the generation of a huge, uni- 
directional current becomes possible. We define the di- 
rected current as the time average of the ensemble aver- 
aged momentum, i.e., 



J 



1 

T 



dt < p(t) > , 



with the ensemble average given by 



1 N 

n=l 



Pn(0))- 



(24) 
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Figure 10: (Color online) The mean escape time vs. the an- 
gular driving frequency is plotted with error bars. The lower 
and upper dashed lines correspond to T m in = T/An — 1/(20) 
and T/4 = it /(2ft), respectively. The number of particles 
emanating from the escape set shown in Fig. [8] is taken as 
N — 512. The remaining parameter values are F = 0.05 and 

e = o. 




Figure 11: (Color online) Temporal behavior of the directed 
current defined in Eq. (124)) for three different angular driving 
frequencies of the inclination-modulation force as labeled cor- 
respondingly on the graphs. The amplitude and the initial 
phase of the external force are set at F — 0.05 and Oo = 0, 
respectively. The number of particles is N = 512. 



Here, N denotes the number of particles constituting the 
ensemble. In Fig. [II] we depict the time evolution of the 
directed current J for moderate, intermediate and very 
slow modulations. The time is scaled according to t — > 
fit; thus assuring that for a given simulation time T s the 
number of oscillation periods of the external force is the 
same regardless of the value of the used angular frequency 
fl. At a moderate driving frequency of fl — 1 the average 
flow (numerically) practically vanishes asymptotically on 
the displayed time scale, whereas for fl = 0.1 the decay 
towards zero proceeds slower. This is in compliance with 
the findings in |l5l |. 



In distinct contrast, however, for extreme slow driv- 
ing at 17 = 0.01, the mean momentum assumes a quasi- 
stationary regime of considerable large size. We note 
that our observed averaged directed current persists over 
a long time interval that agrees well with the explicit es- 
timate for the duration, tf of unidirectional motion given 
in J28(], which adopted to our system notation reads 



*/ 



1 



(\/2^/rj) 5 

W: ln 4 [V2Jr/(Ffl)} 



(26) 



yielding for fl = 10~ 2 and F = 0.05 a figure of the or- 
der of tf ~ 10 8 . Plotting the directed current as as a 
function of the angular driving frequency fl one notices a 
strong decay of the current with increasing angular driv- 
ing frequency (not shown here). The resulting 1/fl de- 
pendence of the current corroborates with the expression 
for the momentum gain given in Eq. 11211) ; see also in Ref. 
[la], [23], providing intuitive arguments that support the 
strong enhancement of the chaotic transport in space. 



VII. CONCLUSION 

In this work we have investigated the dissipation-less, 
time-dependent driven Hamiltonian dynamics of parti- 
cles evolving in a symmetric, spatially periodic potential 
whose inclination is temporally varied periodically by an 
external ac-force. 

We have focused interest on the generation of a di- 
rected flow of ensembles particles which are trapped in 
the interior of a potential well. The choice of the ini- 
tial starting values is subjected to the condition that 
for a static inclination the trapped particles cannot es- 
cape from a potential well. Then, the only possibility left 
for escape from the potential well is due to the chaotic 
dynamics which arises in the system dynamics due to 
the time-dependent forcing term. In fact, upon applying 
the time-periodic modulation of the inclination, trajec- 
tories that become embraced by the developing chaotic 
layer around the broken separatrix may cross the latter. 
However, for fast and intermediate modulation frequency 
there results no substantial directed flow. This is so, be- 
cause the trajectories frequently cross and re-cross the 
separatrix corresponding to leaving and re-entering the 
adjacent wells of the potential wherein the particles dwell. 

For adiabatic inclination modulations we have demon- 
strated that for all initial conditions contained in the es- 
cape set, motion takes place in a unique direction that is 
controlled by the phase of the modulation term. It has 
been shown that the slower the modulation the larger 
is the gain in momentum of the escaped particles and 
thus the emerging asymptotic current that obeys a I/O 
dependence. 

Concerning an explanation of this phenomenon it 
seems that the cantori, partitioning the nonuniform 
chaotic layer, are the less leaky the smaller the modu- 
lation frequency fl. Research regarding the modulation 
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frequency dependence of the sticking times to the bound- 
ary of regular islands is in progress J43| . The cantori form 
almost impenetrable barriers that confine trajectories for 
a very long but transient period. One should remark that 
eventually this transient period of directed motion ter- 
minates because the trajectory escapes through one of 
the holes in the cantori and accesses other regions of the 
chaotic layer. Therefore the motion does not necessarily 
proceed unidirectionally: Unless the trajectory gets cap- 
tured by ballistic channels it itinerates within the chaotic 
layer going along with changes of the direction of motion. 
Suitable physical systems that come to mind to exper- 
imentally test our findings are periodically driven cold 
atom optical lattices, as recently fabricated in study- 
ing dissipative, classical ratchet dynamics [4J, [45|]. In 
order to verify our multi-facetted findings of directed, 
dissipation- less transport the use of off- resonant, far de- 
tuned laser beams is required; thus minimizing the dissi- 



pation in these cold atom set-ups: A scenario proposed 
in different context also for Hamiltonian quantum ratch- 
ets in Refs. [2a, [461]. Moreover, the particle trapping- 
detrapping transitions induced by time-dependent mod- 
ulations of the potential as described in this manuscript 
can also be applied in fluid dynamics to design particle 
traps in incompressible open flows as discussed in [48l |. 

As an interesting extension of the present work we cur- 
rently engage in studying this dissipation-less enhance- 
ment effect for coupled nonlinear systems which are com- 
posed of a chain of interacting particles J47j . 
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